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Abstract 

Light dark matter with mass smaller than about 10 GeV is difficult to probe from direct detection 
experiments. In order to have the correct thermal relic abundance, the mediator of the interaction 
between dark matter and the Standard Model (SM) should also be relatively light, ^ 10^ GeV. If 
such a light mediator couples to charged leptons, it would already be strongly constrained by direct 
searches at colliders. In this work, we consider the scenario of a leptophobic light Z' vector boson 
as the mediator, and study the the prospect of searching for it at the 8 TeV Large Hadron Collider 
(LHC). To improve the reach in the low mass region, we perform a detailed study of the processes 
that the Z' is produced in association with jet, photon, and Z^. We show that in the region 
where the mass of Z' is between 80 and 400 GeV, the constraint from associated production can be 
comparable or even stronger than the known monojet and dijet constraints. Searches in these channels 
can be complementary to the monojet search, in particular if the Z' couplings to quarks {gz') and 
dark matter {go) are different. For gu < gz>, we show that there is a larger region of parameter space 
which has correct thermal relic abundance and a light Z', Mz' ~ 100 GeV. This region, which cannot 
be covered by the mono-jet search, can be covered by the resonance searches described in this paper. 



1 Introduction 



Dark Matter (DM) consists of 24% of the energy density of our universe. However, the nature of it is one of 
the outstanding mysteries. The most popular class of candidates for DM are the stable Weakly Interacting 
Massive Particles (WIMPs), and the annihilation of WIMPs in the early universe determined the observed 
relic abundance of DM. This scenario provides a promising approach to detect the WIMPs directly [1], 
with many recent results [2l410j. The direct detection experiments are more sensitive for WIMPs with 
masses heavier than tens of GeV. Assuming the WIMP-nucleon interaction is spin independent, the limit 
can be as strong as (TwiMP-nucicon < 10^^"^ to 10~^^ cm^. For light WIMPs, the energy transfer during a 
WIMP-nucleus scattering is proportional to the WIMP mass. If the WIMP is sufficiently light, the energy 
transfer is too small to be detected, resulting in a much weaker limit. For example, the upper limit is 
weaker than 10~^^ cm^ for Mwimp ~ 5 GeV. 

Recently, the low mass dark matter scenario has received more attention. For example, there is an 
intriguing connection between the WIMP number density and baryon number density in this scenario [111 - 
I29j . which is motivated by the fact that the energy densities of dark matter and baryons are of the same 
order of magnitude. High energy colliders can play a significant role in the search for light WIMPs. The 
simplest approach is to assume that the interaction between WIMPs and Standard Model (SM) particles 
can be written as effective contact operators |30H37j . with the mediators between the WIMP and the SM 
sector integrated out. The simplest observable signal would be WIMP pair production associated with a 
jet, which is often referred as the monojet + missing transverse energy (MET) search channel. 

On the other hand, it is a generic possibility that the mediator can be light. In this case, it cannot 
be integrated out while considering scattering process at the LHC. Such an effect has been studied in 
Refs. |31H38| . A light mediator is particularly motivated by reproducing the observed thermal relic 
abundance, which requires a WIMP annihilation cross section is of order picobarn. If WIMP annihilation 
mainly proceeds through s-wave , the cross section can be written as 

where we have assumed that the interaction between WIMP and the SM quarks is mediated by a vector 
boson Z'. This is the so called dark portal. The relevant Lagrangian can be written as 

C B gz'qi^qZ'^ + goX^xZ'^, (2) 

where we have assumed that the dark matter particle, denoted by X) is a Dirac fermion. gz' and g^ 
are the coupling of the mediator to quarks and to WIMP, respectively. A'^^ is the number of SM flavors 
that are kinematically available. For such a low mass resonance its coupling to leptons must be strongly 
suppressed, otherwise it would be ruled out by LEP. That's the reason that in our model the Z' couples 
only to quarks in SM (leptophobic). 

Such a light mediator with Mz' ~ 0(100) GeV is within the reach of the Large Hadron Collider (LHC). 
The most straightforward way to search for a leptophobic Z' is a direct search for dijet resonance, a "bump 
hunt", which is summarized in Ref. |38]. Although there is some model dependence, it turns out that 
resonance searches for the vector mediator with masses between 250 GeV and 4 TeV at the colliders and 
mapping to DM direct detection constraint can provide stronger constraints than the monojet searches. 
On the other hand, the constraints become much weaker for low mass mediator. Since the jets from 
lighter Z' decay would be softer, the signal of a light Z' suffers from very low trigger efficiencies. 

In this work, we study a class of alternative channels in which the Z' boson is produced in association 
with a hard jet, a hard photon, or a massive gauge boson in the SM. We also consider the case of Z' pair 
production. The rates of the SM background for the associated productions can be significantly lower. 
Therefore, we can lower the trigger threshold and enhance the signal efficiency. In the following we will 
demonstrate that this is indeed a promising way of finding or constraining the leptophobic Z' mediator, 
and eventually the DM particle. 
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We present a detailed description of the Monte Carlo simulation in Sec. 2. In Sec. 3, we scan the 
parameter space and show the reach for each channel and their combination. In Sec. 4, we relate the Z' 
search to the collider search of dark matter. We summarize the results in the Sec. 5. 



2 Simulation and Event Selection 
2.1 Signal and Background Simulation 

The signal and background are both generated by MadGraph/MadEvents 5.1.3.9 [39]. We used the kT 
MLM Matrix Element-Patron Shower (ME-PS) matching algorithm. The decay width of Z' is calculated 
by BRIDGE 2.24 [50]. The following decay and hadronization procedure are performed by pythia-pgs 
2.1.8. The jets are constructed by f ast j et 3.0.0 [51] with i? = 0.6 anti-Zcy jet algorithm, and smearing 
is performed for jets, leptons, photons with the CMS-like energy resolution 

bEj 0.84 



0.07, (3) 



6Ei 0.028 0.125 

-ET- = I , ® 0-003 ^ (4 
El ^/E^/GeV Ej/GeY' 

0.007. (5) 



6E^ 0.01 



E^ ^E^/GeV 

Throughout this paper we assume a 15 fb^^ total integrated luminosity at centre-of-mass energy of 8 
TeV. 

The next-to-leading order correction K-factor K = (Tnlo/(^lo will not significantly change our results. 
In particular, it is calculated by MCFM 6.2 [53] to be roughly one, both for the Z' + channel and the 
Z' + Z^ channel. For the Z' + jet and the Z' + 7 channels we do not preform a calculation, but |44j 
suggests the K factor for pure QCD 3 jets should be about 1.3 or so, and in [55] the K factor for 7+ 2 jets 
is about 1.2 or so. All the K factors are close to one, so we expect counting the higher order corrections 
will not significantly change our reach. 

We approximate Z' as a Breit-Wigner resonance. The decay width of Z' can be written as 

T^^Mz'{NfNc + r^), (6) 
Izvr 

where Nj = 5 or 6 if Mz' is above the ditop threshold, Nc = 3 is the number of colors in QCD, and 
r = Qd/qz'- As we will see in the following sections, the narrow width condition Tz' ^ Mz' is always 
satisfied in the parameter region we are interested in. Therefore, the Breit-Wigner approximation is valid 
in our case. We will not consider the interference effect of the Z' signal with the SM background. Such 
an effect is only important within the Z' width. So in our narrow width Z' scenario the effect is small, 
and it will be completely washed out due to the jet energy resolution. We verified that this is indeed the 
case. 



2.2 Event Selection 

The most important class of event selection cuts for our signal are the jet energy thresholds and acceptance 
cuts. We begin with the Z' + jet channel. Due to the large QCD background, we have to adopt relatively 
high jet energy threshold. With our choice, the leading order SM background cross section corresponds to 
an event rate of a few Hz or so, with the current peak instantaneous luminosity. Since the energy of the 
jets from Z' decay is closely correlated with the Z' mass, we use two complementary sets of pTj selection 
cuts, so that in combination they give us good signal efficiency for a large range of Mz'- For light Z' 
{Mz' < 350 GeV) we require 

PTji > 350 GeV PTJ2,PTJ3 > 70 GeV. (7) 
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Z' Mass (GeV) 


60 


80 


100 


120 


150 


250 


350 


Mjj (GeV) 


50-70 


60-90 


80-110 


100-130 


120-160 


210-270 


290-370 


Z' Mass (GeV) 


450 


550 


650 


750 


850 


950 




Mjj (GeV) 


370-480 


450-590 


530-700 


610-800 


700-900 


800-1000 





Table 1: Invariant mass window for each hypothetical Z' mass. We optimize this with signal Monte Carlo, 
where the boundaries correspond to the bins with half of the peak bin height. 



In this case, the second and third hardest jets are dominantly coming from Z' decay. We choose to use a 
relatively low threshold for them to enhance signal efficiency. At the same time, we have to require the 
hardest jet to be very energetic to suppress the rate of QCD background. For heavier Z' [Mz' > 350 
GeV), the decay products of Z' give two hard jets that are most likely to be the two leading jets. The 
two hardest jets, if coming from Z' decay, are typically close in py. Therefore, equal thresholds for the 
first two leading jets are desirable. We also require the jet produced in association with Z' to be energetic 
to suppress the QCD background. In particular, we require 

Ptj > 160 GeV for three leading jets. (8) 

In the other channels, Z' -|- 7 and Z' -|- jZ'^ ^ the signal also contains at least two jets coming from 
Z' decay. At the same time, the additional hard object in the event, such as a hard photon or hard leptons 
from jZ'^ decay, can be used to efficiently trigger on this class of events. Therefore, we can afford to 
use lower thresholds for jets 

Vtj > 50 GeV for two leading jets. (9) 

We will focus on central jets with a good energy resolution. For all the channels considered in this 
paper, we impose an acceptance cut 

hj\ < 2.5. (10) 

We apply this cut not only to jets but also to 7s and charged leptons. 

We also note that the two jets from Z' decay tend to be central. At the same time, the background 
QCD jets are more forward with larger rapidity gaps. Therefore, we impose an rj separation cut 

|Ar/jj| < 1.7, (11) 

for the two jets which are used to reconstruct the Mz' ■ 

We have also considered different jet energy thresholds in both light and heavier Z' cases. The results 
are tabulated in Table H] and Table [5] in Appendix \^ While we use the relatively conservative choices of 
the jet energy threshold, lower them can certainly increase the reach. 

To search for the Z' resonance, we have to identify two jets as candidates for Z' decay products. In 
the Z'+^et channel with Mz' < 350 GeV, we choose the second and third hardest jet. For A'lz' > 350 
GeV, we choose the two leading jet instead. In Z' + j and Z' + W^/Z^ channels, we also choose the two 
leading jets. After identifying the two candidate jets, we require their invariant mass to be within the 
mass window around the target Z' mass. The width of the window has two origins, one is the natural line 
width or the physical Z' total decay width, and the other is due to the detector finite energy resolution. 
In our case, the latter one always dominates. We optimized the mass window for a set of Z' masses 
based on Monte Carlo simulation. The complete list of the invariant mass windows used in our analysis 
is Table [221 

We also impose selection cuts on the other (non-jet) hard objects in the signal. For the Z' + ^ channel, 
we require 

PT-y > 50 GeV. (12) 
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For the Z' + channel, we focus on the leptonic channel. We require 

PT^ > 25 GeV, (13) 
> 25 GeV. (14) 

For the Z' + Z^ channel, we consider two different Z^ decay modes. The first channel is the charged 
leptonic Z^ decay mode. In this case, there are two opposite sign electrons or muons which correctly 
reconstruct a Z^. Specifically, we require 

Pti > 25 GeV for two opposite sign leptons, (15) 
85 GeV < Mu < 97 GeV. (16) 

The second channel is the invisible decay Z^ > uv. We require 

> 60 GeV, (17) 

and veto any charged leptons. They should be replaced by the practical triggering conditions at the LHC. 

In our study, we have not considered Z' — )■ 66 as a possible decay channel, mainly for simplicity. The 
results are already quite encouraging without 66 decay channel, and including it will certainly enhance 
the discovery reach. Extending the analysis in this paper to this case is straightforward, after properly 
taking into account the b-tagging. For heavier Z', decaying into would also give another signal channel. 
At the same time, identifying the top requires different strategies, depending on Mz' 06]. 



3 Reach of Different Channels 

In general, the couplings of Z' to the left and right-handed quarks can be different. In the Z' + jet and 
Z' + 7 channels, since QCD and QED are vector-like, only the combination g = {g^/L + 9z'_r)^''^ 
relevant, where gz'L and gz'R are the couplings of Z' to the left and right-handed quarks. Only the left 
handed coupling is relevant for the Z' + channel. The most complicated channel is Z' -\- Z^, in which 
both the gz'L and gz'R are relevant. Further complication comes with two kinds of decay channels under 
consideration. However, we will see in the following that the reach is mainly due to the non-chiral Z' + 
jet channel. The chiral Z' + contribution is nearly the same as the non-chiral Z' + j contribution. 
Z' + Z^ has a very small contribution. Therefore, we expect the ratio gz'i/dz'R will not play a significant 
role, and we will only show the results of gz'L = gz'R = dz'- 

As illustration, a list of cross sections for various background channels and signal channels with 
gz' = 0.35 are shown separately in Table [2] and Table [3l The expected 2a constraints from different 
channels are shown in Fig. [U where the red, green, blue, purple and black curves correspond to the upper 
limits from the Z' -|- jet, Z' + j, Z' + W^, Z' + Z^ channels and the combined constraint of the four 
channels, respectively. Here we have assumed here go = 0. In realistic models considered later, the 
decay branching ratio to jets will induce a suppression to the signal rate in each channel, however, the 
suppression is small at least for gj^ ^ A detailed discussion is in Sec. [H For Mz' > 1000 GeV the 

background QCD dijet rate is low enough that the more efficient direct dijet resonance search provides a 
better limit [59], so we only consider the cases with Mz' < 1000 GeV here. 

From Fig. [T] we can see that the constraint from Z'-|-jet channel is stronger than other channels, 
especially in the large Mz' region. The dominant reason is the QCD coupling is much larger than the 
other relevant gauge couplings. For the Z' + ^ channel, both the signal and the background are suppressed 
by the fine structure constant of electromagnetic interaction, aem- Whereas in the case of Z'-l-jet channel, 
it is replaced by as, which is much larger. As a result, both the signal and the background increase by 

^Note that due to the large number of degrees of freedom in the quark sector, and our assumption that there is only one 
species of dark matter, go = gz' only leads to a small invisible width, BRinv ~ 1/16. 
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QCD 3 jet 


dijet+7 


dijet+TF^^ 


dijet+(z" ^ e+r) 


dijet+(Z" vv) 


PTJ, V cuts 


487 pb 
412 pb 


327 pb 


160 pb 


7.6 pb 


20 pb 


At] cut 


370 pb 
284 pb 


zzD pb 


iuy pb 


4.Z pb 


14 pb 


Mjj £ [120 - 160] GeV 
Mjj E [290 - 370] GeV 
Mjj G [800 - 1000] GeV 


42 pb 
67 pb/11 pb 
25 pb 


49 pb 
18 pb 
0.76 pb 


13 pb 
4.0 pb 
0.21 pb 


0.95 pb 
0.32 pb 
0.04 pb 


2.5 pb 
1.3 pb 
0.11 pb 



Table 2: Leading order standard Model background rate for each of the channels, after various selection 
cuts. In the first row the basic pT thresholds and rj < 2.5 cut for all particles in an event are implemented. 
Then in the second row we add Ar]i2 cuts for two leading jets, as well as the invariant mass window for 
e^e~ or jjL'^ in the dijet+ {Z^ — )■ channel to reconstruct . The event selection cuts are 

designed in particular to make the corresponding rates here no more than a few Hz for each channel. 
Lastly invariant mass windows are applied for three hypothetical Z' masses. In the QCD 3 jet channel, 
we list the two cross sections, which correspond separately to the 350-70-70 pTj triggering (the former) 
and the 160-160-160 pTj triggering (the latter). 





Z'+ jet 


Z' + -i 


Z' + 


z' + {z'^ ^ i+r) 


Z' + {Z^ u9) 


Mz' = 150 GeV 
Mz' = 350 GeV 
Mz' = 950 GeV 


0.80 pb 
0.72 pb/0.15 pb 
0.20 pb 


0.76 pb 
0.30 pb 
0.013 pb 


0.45 pb 
0.14 pb 
0.006 pb 


0.031 pb 
0.0096 pb 
0.0004 pb 


0.054 pb 
0.021 pb 
0.0011 pb 



Table 3: Leading order Z' signal rate for each of the channels. All selection cuts in Table 2 (including 
the mass window) are applied. Here we are assuming Ql = Or = 0.35. Again in the Z' +]et channel we 
list the two cross sections, which correspond separately to the 350-70-70 pxj threshold (the former) and 
the 160-160-160 pTj threshold (the latter). 



roughly the same amount in comparison with the Z' + 7 channel, and therefore the S/^/B for the Z'+jei 
channel is larger than that for the + 7 channel. This is exactly what happens in the large Mz' in 
Fig. [TJ In the small Mz' region, as shown in Tables [2] and [3l the cuts make the cross sections of Z'+jei 
channel and Z' + 7 channels comparable in both signal and background processes. Therefore, as shown 
in Fig. [H the constraints on gz' is comparable in the region Mz' ~ 150 GeV. 

In the Z'+jet channel, in the small Mz' region where Mz' ^ 100 GeV, to balance the transverse 
momentum of the hardest jet, the Z' resonance is highly boosted. Therefore, the two jets from the decay 
of the Z' resonance are close to each other, and cannot be distinguished from a single fat jet. Due to the 
same reason, the small M23 region in the background is also removed, where M23 is the invariant mass 
of the second and third hardest jets. As a result, the events of both the signal and the background are 
cut off by the requirements in Eq. ([Tj) and the statistics becomes poor, and therefore the limit on gz' 
from Z'+jet channel weakens significantly in the region where M^/ < 80 GeV. A quantitative analysis 
which shows the scale of this effect can be found in Appendix [Bl For example, for M'^ ~ 80 GeV which 
is discussed in Appendix [Bl if we use the R = 0.1 anti-fer jet algorithm rather than the default R = 0.6 
one in a hadron level analysis, we will get about 20 times more events with a dijet invariant mass in the 
window of 60 — 90 GeV. Jet substructure techniques should be able to help in this regime, and we leave 
the details in this direction for a future study. In contrast, for the Z' + 7 and Z' + W channels, there is 
no such kinematic configuration. Therefore, in the small M^/ region, the constraint is mainly from these 
two channels. 

In the Z' + 7 channel, it is possible for jets to fake photons. For example, there can be a hard jet 
which accidentally becomes a vr*^ and passes the neutral electromagnetic trigger to mimic a photon. The 
jet fake gamma rate is tiny (less than 0.3%), but given the QCD jets are much more abundant, the fake 
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Figure 1: Reach plots for S/^/B for the Z' + jet (red), the + 7 (green), the Z' + (blue) and 
the Z' + Z^ (purple) channels, for 15 fb~^ integrated luminosity and S/y/B = 2 which corresponds to 
95% confidence level. All the channel is assuming gz'L = 9z'R = Qz', and in the Z' + channel the 
reach from two charged leptons decay product and from two neutrinos are combined. Eventually their 
combination for exclusion is shown as the black thick line. 

trigger is still considerable. Taking this into consideration, the background will increase by about 20%. 

The Z' + channel has comparable S/^/B ratio with the Z' + 7 channel over all the Z' mass 
region, which is just a coincidence of various physical reasons contributing in different directions. For 
example, the weak coupling is larger than the electromagnetic coupling, the is massive so that the 
production is relatively suppressed. — )• £^ has further suppression from the leptonic decay branching 
ratio. At the same time, the 7s are concentrated in the collinear region and hard to pass the pT^y > 50 
GeV selection cut, while the acceptance of a — )• £^ is higher. 

For the Z' + Z*^ channel, we have to combine its two decay channels. The neutrino channel has a 
larger branching ratio (~20%) than the chagerd lepton channel (~6.7%), and the former has a slightly 
better S/^/B. However, even the combined signal significance is much smaller than other channels, so 
this channel is less interesting. 

We have also checked the Z' + Z' pair production channel. The S/^/B ratio is always much less than 
1 in the region of mass and coupling we focus on. One reason is that the couplings are all relatively small, 
which leads to small production cross section. At the same time, the signal is in a pure 4-jet final state, 
which is overwhelmed by the dominant QCD background. For a similar background rate with the 3-jet 
case we have to use nearly the same jet pTj thresholds, but the signal cross section is further suppressed 
by a small factor of g^^,. Therefore, the reach in this channel would be much weaker, and we will not 
provide the full analysis here. 

Fig. [2] shows the comparison between the constraints from dijet search with associated products, 
monojet search and direct dijet searches. To make connection with dark matter detection, we will now 
consider go ^ Q and introduce r = go/gz' to parameterize its size. As an illustration, we assume 
r = go/gz' = 1 in Fig. [2]. The red curve is the 95% C.L. upper limit from ATLAS monojet searches |47j . 
the green and blue curves are the 95% C.L. upper limit from the CDF [38] and ATLAS [l9] dijet resonance 
searches, and the purple curve shows the 90%C.L. upper limit from UA2 dijet resonance search |50j . 
respectively. For dijet searches, colliders with smaller centre-of-mass energy give stronger constraints, 
since when Mzi is much smaller than the centre-of-mass energy of the collider, the constraint suffers 
from large QCD background due to the peak of the gluon parton distribution function at low x. For the 
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Figure 2: Comparison between the combined constraint of different channels shown in Fig. [T] and the 
constraints from dijet and monojet searches. We have assumed r = go/g'z — 1- The thick black curves 
shows the combined constraint as in Fig. [TJ The red curve corresponds to the 95% C.L. ATLAS monojet 
upper limit. The green and blue curves correspond to 95% C.L. upper limits from CDF and ATLAS dijet 
searches, respectively. The purple curve corresponds to 90% C.L. upper limits from UA2 dijet search. 



comparison between the associated dijet constraint and the monojet constraint, a large amount of the 
background events can be removed with the help of the invariant mass window cut, so the associated 
dijet constraint can be stronger than the monojet constraint. On the other hand, as discussed before, 
in the very light Z' region {M'^ < 80 GeV), the two jets from the decay of Z' are either highly boosted 
and cannot be distinguished from a single jet, or probably cut by the pxj > 50 GeV threshold which is 
roughly half of the resonance mass, whereas the invisible decay of Z' is only characterized by large missing 
transverse energy. Therefore, in this region the monojet constraints can be stronger than the associated 
dijet constraint. 

Before the end of this section, we briefly mention some existing results in the + dijet resonance 
channel. Mainly motivated by checking the CDF + dijet anomaly, the ATLAS [51] and more recently 
CMS group [52] have performed searches in the same channel. In [51] based on 5 fb~^ integrated luminosity 
a leptophobic Z' of 150 GeV and gz' — 0.2 [53] is excluded. This is in broad agreement with our results 
shown in Fig. [TJ Possible strategies to enhance the LHC reach in the + dijet has also been studied 
in [53]. However, the kinematics of this model in which the djiet resonance and W are decay products of 
a heavier new resonance, is very different from the scenario considered in this paper. 



4 Z' as a. Portal between the SM and Dark Matter 

The Z' can mediate interaction between dark matter and SM particles, forming the so called dark portal. 
In this case, the constraint on gz' can be mapped onto the constraints on DM direct detection cross 
section. The direct detection cross section for a nucleon (proton or neutron) is 

9gl,glMj,Ml ^ ^g_z^Y .^.2 / lOOGeV y , 



^si - -^f^^J-N-u ^ 7.7 X 10-^" ^ cm^ (18) 

where M^v is the mass of the nucleons, and M/j = 5 GeV is assumed. 
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Mz(GeV) 

Figure 3: Comparison of the gz' = go reach of our Z' association production (soUd and dashed black 
curves for 2a expected constraints and 5(T reach, respectively) with other experiments. The red curve is 
the constraint from ATLAS monojet search with 10 fb~^ integrated luminosity at 8 TeV. The green curve 
is the bound from dijet resonance search by CDF with 1.13 fb~^ integrated luminosity. The blue curve is 
the constraint from ATLAS dijet resonance search with 1 fb~^ data set. 



For gz' = go, the constraints on direct detections cross section are shown in Fig. [3l The major 
improvement is in the region with Z' lighter than limit from the CDF dijet pole search. The constraint 
can be as strong as a few xlO~^^ cm^. Assuming gz' = go and Mz' > 80 GeV, limits from associated 
production are also stronger than those from the ATLAS monojet search. The current bound assumes an 
integrated luminosity of 15 fb~^, this constraint will becomes stronger if the assumed integrated luminosity 
increases. 

Relaxing gz' = go leads to interesting scaling behavior. The production rate of Z' is proportional 
to g^,. The decay branching ratio of Z' into dijet final states in the case of gz' 7^ go can be written as 
g'^iN fNc / {g'ziN fNc + gjj). On the other hand, the decay branching ratio of Z' into DM can be written 
as g'jj/{g'^,NfNc + gjj)- Therefore, for a general r = gn/g'z^ the cross sections for monojet and dijet 
processes can be written as 

(0) Nc'Nf + 1 2 NcNf + 1 

^monojet (r) = a„j,, X ^^^^^ oc gz'gD X r ; 



NcNf + " " NcNf + 

(0) NcNf + 1 1 NcNf + 1 

^ ^ oc qz'QD X rr 

NcNf + r2 r NcNf + 



where Cmonojot ^^'^ "^dijct cross sections for r = 1. Therefore, as long as the Z' is narrow-widthed 

and light enough so that it can be produced on shell, one can get the constraints on the coupling for a 
general value of r by scaling the constraints shown in Fig. [2] using Eq. (jl9p . 

For fixed gz'dD, the constraints of the production cross section from monojet and dijet searches have 
approximately opposite dependence on r. Therefore, the monojet method and our associated dijet method 
are complementary: the former works for large r whereas the latter works for small r. We can also take 
into account the theoretical consideration of DM thermal production, or the so called "WIMP miracle". 
If the Z' is the only portal between the DM and the SM sectors, and the relic abundance of DM is 
determined by thermal freeze-out, the DM annihilation cross section to SM particles through the Z' will 
be given in Eq. ([1]), where we can see that (av) depends on the couplings only through the product gz'9D 
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Figure 4: Comparison of the 5a reach for the discovery of a 5 GeV WIMP by the Z' associated production 
search (all solid curves) and monojet search (all dash-dotted curves), for different ratios of r = go/dz'- 
The region of parameter space above a certain curve either has been excluded or can be constrained 
by the corresponding search channel. All curves in red correspond to go/dz' = 0.3, all curves in green 
correspond to gn/dz' = 0.7, and all curves in blue correspond to go/gz' = 1- In shaded region, the 
XX ^ process generates the correct relic abundance, namely requiring the cold dark matter density of 
the universe Vt^dmh'^ = 0.111 ± 0.006 [55]. 



in the case that Md <^ Mz'- To get the observed relic abundance, {av) is fixed to be around 3 x 10~^^ 
cm'^s"^; as a result, the product gz'go becomes a function of Mz' for fixed M^, which is shown as the 
grey band in Fig. [H 

In the region that Mz' and r are both small, the monojet constraint becomes weaker, as expected 
from Eq. ()19p . It cannot probe the shaded region in Fig. [4] where the correct relic abundance is generated. 
In the region r < 0.7, from the crossing of the solid and dash dotted green curves in Fig. [U we can see 
that the correct relic abundance can be generated with parameters that satisfy the monojet constraint 
and can be reached by the associated dijet search with a significance larger than 5a at the same time. In 
the region of Mz' < 130 GeV, from the red curves in Fig. U we can see that if r is smaller than about 
0.3, the associated dijet resonance search can reach the full region of parameter space with correct relic 
abundance, with more than 5a significance. Whereas in the region of Mz' > 130 GeV, the associated 
dijet bump search cannot reach 5a significance due to the constraint from UA2 dijet resonance search as 
shown in Fig. [2j 



5 Summary and discussions 

In this paper, we have suggested a new systematic way of finding/constraining the general leptophobic 
Z' gauge boson. To overcome the large QCD background at the LHC, we considered the channels in 
which Z's are produced in association with a hard jet, an energetic photon, or a W^/Z^. We performed 
a detailed study of the potential of LHC search in these channels, and show that the Z' + jet, + 7 
and the Z' + channels are independently promising, and a combination can significantly improve the 
reach. 

Assuming this Z' also couples to dark matter, forming the so-called dark portal, this bound can be 
mapped onto a dark matter direct detection bound. We demonstrate that there is improvement in the 
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low Z' mass region (M^/ < 400 GeV) from the use of associated production. The cohider constraint of 
DM direct detection cross section can be extended to a smaher hypothetical Z' mass. For gz' > Qd^ the 
Z' associated production method holds an advantage over the monojet method. For example, if the Z' 
mediates the only channel for DM to annihilate into SM particles, with small value of r, the monojet 
constraints become weak and the parameter space for generating the correct relic abundance with small 

5 GeV Md and weak scale Z' is still allowed, as shown in Fig. [H At the same time, large region of this 
parameter space can be covered by the associated dijet bump search discussed in this paper. In particular, 
we can discover such a Z' in this search channel with a significance larger than 5fJ, if go/dZ' is smaller 
than about 0.7. 

6 Acknowledgments 

We would like to acknowledge useful information from Antonio Boveia about the triggering of the ATLAS 
experiment. R.H. wish to thank Peter Skands and Xiaohui Liu for useful discussions. H.A. is supported by 
in part by the Government of Canada through NSERC and by the Province of Ontario through MEDT. 
L.T.W. is supported by the NSF under grant PHY-0756966 and the DOE Early Career Award under 
grant ,de-sc0003930. 

A Variations on jet energy thresholds 



Threshold 


350-70-70 


300-70-70 


400-70-70 


350-60-60 


350-80-80 


Background 


370 pb 


767 pb 


200 pb 


410 pb 


333 pb 




15.0 


19.4 


8.9 


15.2 


14.8 



Table 4: Comparison of small variation of the first set of jet energy thresholds for the 3-jet signal. For 
the leading order event rate we also impose the |A?723| < 1-7 cut, but no invariant mass window choice. 
In S/^/B we are assuming gi = gn = 0.35, Mz' = 150 GeV and 15 fb~^ integrated luminosity. 



Trigger level 


160-160-160 


150-150-150 


140-140-140 


170-170-170 


180-180-180 


Background 


286 pb 


406 pb 


579 pb 


205 pb 


148 pb 


s/Vb 


10.8 


11.2 


12.6 


8.9 


8.0 



Table 5: Comparison of small variation of the second set of jet energy thresholds for the 3-jet signal. For 
the leading order event rate we also impose the |A?7i2| < 1.7 cut, but no invariant mass window choice. 
In S/^/B we are assuming gz = SR = 0.35, Mz' = 450 GeV and 15 fb~^ integrated luminosity. 

We compare different choices of jet energy thresholds, and the results are tabulated here. 

B Low Mz' detection thresholds for different channels 

We start with considering an event with three jets in the final state, labelled by Ji, J2 and J3, with 
transverse momenta pri, Pt2 and pTs, respectively. Here, we require that pxi > Pt2 > PT3- Then, the 
invariant mass of J2 and J3 can be written as 

Mis = 2(p|.2 + p'2)'/'(PT3 + P's)'^' - 2PT2 • PT3 " 2p2.P3. , (20) 

where we neglect the invariant mass of a single jet. p2z and p^z are the longitudinal momenta of J2 and 
J3, respectively. We can always work in a frame where Pz3 = such that the above expression can be 
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simplified as 



^23 = '^iPT2 + P^2)^^^PT3 " 2pT2 • PT3 

Pz2 \ 



2pT2PT3 



1+ 2 / 

Pt2J 



+ 2p2TP3t(1 - COS (/)) , (21) 



where (/> is the angle between andpsr- It is easy to see that 2pT2PT3 > 2p™3°(p™"^— p^g'^) ^ (200GeV)^, 
where p™}", and p™3° are the cuts we impose on pxi , Pt2 and pxs ■ In practice, for the region AIz' < 350 
GeV, we choose pf^ = 350 GeV, p^^ = pff = 70 GeV. Therefore, in the region Mfg < (200GeV)2 we 
have 

Ml « 2PT2PT3 f 1^ + ^</'' 
V ^Pt2 ^ 

~ PT2PT3(y' + , (22) 

where y is the rapidity of J3. In the anti-ky algorithm, the distance between J 2 and J3 is defined as 

d23 = min (p-„p-3) « -3-—^ . (23) 

If c?23 < 1/Pr2' and J3 will be identified as a single jet. Therefore, it requires ^23 > l/pj'2 *° identify 
J2 and J3 as two separate jets, which means 



M23 > RVp^Fp^ > R^P^TWi - PtT) - 80 GeV , (24) 

where R = 0.6 has been used. Therefore, the cuts we imposed on the three-jet final state configuration 
remove both the background and the signal in small M23 region as shown in Fig. [5l We should note 
that this argument is based on parton level analysis. Taking the parton shower and realistic detector 
effects in to account, a small amount of events from both background and signal in the region M23 < 



R\l PtT^PtT ~ Pt™) niay pass the cuts and leak into the signal region. However, the statistics in the 
small M23 region remains very poor, and so does the constraint on the coupling. 

For the Z' + 7 channel, as shown in Fig. [5l since we don't require a very large transverse momentum 
for the photon, the limit of the invariant mass of the two jets can simply be written as 

M12 > 2p^JR/2 ^ 30 GeV , (25) 

where p™™ = 50 GeV is used. Therefore, the constraint from this channel dominates over the Z'+jet 
channel in the small Mz' region. 

For the Z' + W channel, since we only consider the leptonic channel, W boson does not have to be 
produced with a large transverse boost to pass the selection cuts. Therefore, Z' don't have to be very 
boosted, and there will be no limitation on M12 in the small Mz' region. For the Z' + Z'^ channel, however, 
the dominant decay channel of Z^ we are considering is to a pair of neutrinos and the > 60 is required, 
which indicates that the pT of boson must be larger than 60 GeV. For fixed dijet invariant mass M12, 

1 /2 1 /2 

in the Z' + W channel, sj^^^ = Mw + M12, where sj^^^ is the minimal parton level centre-of-mass energy 
to make this process possible. Whereas in the Z' + Z° channel, s^j^ = \J M| + p™™^ + \/ M"^, +p^^z'^. 

1 /2 

Therefore, for both signal and background, in the small A'lz' region, sj^^^ decreases faster in the Z' + W 
channel than in the Z' + Z^ channel with the decreasing of Mz'- As a result, the parton luminosity 
increases faster in the Z' + W channel than in the Z' + Z^ channel with the decreasing of M^/. Thus, 
the upper limit on gz' from the Z + Z' channel increases more rapidly in the small Mz' region with the 
decreasing of Mz' as shown in Fig. [TJ 
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Dijet invariant Mass (GeV) 

Figure 5: Background distributions of dijet invariant mass in Z'+jet, Z' + 7 and Z' -'rW channels. For 
Z'+jet channel, it is the invariant mass of the second and third hardest jets shown in the plot. 
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